Test of CPT and Lorentz Invariance from Muonium Spectroscopy 
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Following a suggestion of Kostelecky et al. we have evaluated a test of CPT and Lorentz invariance 
from the microwave spectroscopy of muonium. Hamiltonian terms beyond the standard model 
violating CPT and Lorentz invariance would contribute frequency shifts &u\2 and 5v2,4 to v\2 and 
1*34, the two transitions involving muon spin flip, which were precisely measured in ground state 
muonium in a strong magnetic fleld of 1.7 T. The shifts would be indicated by anti-correlated 
oscillations in u\2 and v-n at the earth's sidereal frequency. No time dependence was found in vi2 
or z/34 at the level of 20 Ifz, limiting the size of some CPT and Lorentz violating parameters at the 
level of 2 X 10"'^^ GeV, representing Planck scale sensitivity and an order of magnitude improvement 
in sensitivity over previous limits for the muon. 
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Much current theoretical work is devoted to finding a 
more fundamental and general underlying theory from 
which the standard model of particle physics could be 
deduced as the low energy limit. String theory is a cen- 
tral candidate, with the feature that the assumptions un- 
derlying the CPT theorem are invalid for strings, which 
are extended objects. CPT violation then becomes a pos- 
sible signature of more fundamental underlying theories, 
which can be probed by sensitive experimental tests of 
CPT invariance. 

Some years ago Kostelecky and coworkers [1-3] devel- 
oped a plausible extension of the standard model based 
on spontaneous breaking of Lorentz and CPT symme- 
try in an underlying theory without gravity. Their low 
energy effective theory provides a theoretical basis for es- 
tablishing quantitative bounds on CPT invariance. The 
analysis was done in the context of conventional rela- 
tivistic quantum mechanics and quantum field theory 
in four dimensions, retaining the usual gauge structure 
and renormalizability. The Lorentz and CPT violating 
additions to the standard model Lagrangian are highly 
suppressed to remain compatible with established exper- 
imental bounds. 

As applied to muonium (/x^e^ bound state), the most 
relevant terms in the extension to the standard model 
are those in the QED limit, involving only muons, elec- 
trons and photons. The additional terms in the Lorentz- 
violating Lagrangian lead to a modified Dirac equation 
and are given by [2] 
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violating, while a and b are CPT odd, and H, c and 
d are CPT even. The lepton number violating terms al- 
lowed in the extension are not relevant here and have 
been omitted. 

To predict the perturbations to muonium, we consider 
the energy level diagram of the ground state of muonium 
in a magnetic field (Fig. 1). The transition frequencies 
1^12 and z^34 have been measured [4] with high precision 
and are used to determine the hyperfinc structure interval 
Av and the ratio of the muon magnetic moment 

to the proton magnetic moment. 

Leading-order Lorentz violating energy shifts 5vi2 and 
(5j^34 due to the new terms in the Lagrangian can be cal- 
culated using perturbation theory and relativistic two- 
fermion techniques. For our observed transitions at the 
strong magnetic field of 1.7 T, dominantly only muon 
spin flip occurs so these energy shifts are characterized 
by the muon parameters alone of the extended theory. 
This approach results in : 



5vi2 



5vzA w 63 /tt, 
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where 63 — feg 
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The lepton fields are denoted by ^ = e , , and iDa = 
ida — qAa where q = —\e\. All terms are Lorentz- 



^ -f- -fffj [5] are laboratory 
frame parameters. High precision experiments on muo- 
nium (Af) can measure or set limits on the parameters 
of these symmetry violating terms, which are sensitive at 
the Planck scale level [5] . 

Predictions of the values of 1^12 and v^i from standard 
theory - dominantly the QED terms - requires values for 
many atomic constants including m^, /i^, and Av as well 
as the calculation of higher order QED radiative correc- 
tions. The relevant constants and calculations are not 
known to as high accuracy as the experimental determi- 
nations of i>i2 and 1^34. Indeed several of these constants 
- Al/ and lip,/ lip - are obtained from the muonium ex- 
periment. Comparing predictions for V12 and 1^34 (based 
on independent determinations of the required atomic 
constants) with the experimental results results in poor 
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sensitivity to the non-standard model energy shifts ^1^12 
and Si'34. 

However, the theory with CPT and Lorentz violation 

involves spatial components in a celestial frame of refer- 
ence, and since the laboratory rotates with the earth, 
these spatial components vary with time, and conse- 
quently the experimentally observed i>i2 and may 
oscillate about a mean value at the earth's sidereal fre- 
quency O = 27r/23 hr 56 m with amplitudes Si'12 and 
(5z^34. No such signal would be obtained from the stan- 
dard model. In the non-rotating celestial frame of refer- 
ence with equatorial axes {X, Y, Z} where Z is oriented 
along the earth's rotational North Pole, an experimental 
constraint on 6^12 implies [5] 

i|sinx|V(^^)' + (&^)' < (5) 

in which x ~ 90° is the angle between Z and the quanti- 
zation axis defined by the laboratory magnetic field at 
Los Alamos where the muonium experiment was per- 
formed. The transformation from the lab frame quan- 
tity 63 to the non-rotating celestial frame quantities bj 
(where J = X,Y, Z) is given by 

6^ = 6^cosx + {b'^ cos fit + b'^r sin fit) sin X- (6) 

The non-rotating frame quantities, bj are defined by bj = 
bj + m^djQ + ^^jklH'^]^ [5]. The experiment has poor 
sensitivity to the celestial frame parameter 5^, but ideal 
sensitivity to the pair 6^ and by- 

The sum of the transition frequencies, 1^12 -|- 2^34, is 
equal to the ground state hyperfine splitting, Ai/, and 
since we expect (see Eqn.4) 61/12 + Suu ~ 0, no sidereal 
variation is expected in the hyperfine interval. At the 
high field strengths of this experiment, the difference in 
transition frequencies 1^12 — ^34 is almost proportional to 
the magnetic moment of the muon. By looking for a 
sidereal variation in the difference 61^12 — di'34 ~ 263/71, 
we are essentially probing a possible sidereal variation of 
the magnetic moment of the muon (while the g factor 
remains constant). 

The accurate measurements of 1/12 and 1/34 were done 
in a microwave magnetic resonance experiment [4] . A4uo- 
nium was formed by electron capture by muons stopping 
in a krypton gas target. Resonance lines were observed 
by varying the magnetic field with fixed microwave fre- 
quency and by varying the microwave frequency with 
fixed magnetic field. A line narrowing technique was 
used involving observation of a transition signal only 
from M atoms which have lived considerably longer than 
~ 2.2 /xs (Fig. 2). The values reported for 1/12 and 
z/34 at a magnetic field strength corresponding to a free 
proton precession frequency of 72.320 000 MHz were 

i/i2(exp) = 1 897 539 800 (35) Hz (18 ppb), (7) 
2^34 (exp) = 2 565 762 965 (43) Hz (17 ppb) (8) 



in which one standard deviation errors include both sta- 
tistical and systematic errors. To search for a time de- 
pendence of 1/12 and 1/34, we have employed the following 
algorithm. Data from each resonance line run (each last- 
ing about half an hour) are fit at the measured magnetic 
field strength and Kr pressure to determine provisional 
line centers for z/12 and 7/34. These line centers were then 
transformed to their values in a magnetic field strength 
corresponding to a free proton precession frequency of 
72.320 000 MHz. The data were taken at Kr pressures 
of 0.8 and 1.5 atm, so the line centers were corrected for 
a small quadratic pressure shift, and then were extrap- 
olated linearly to their values at zero pressure, using a 
pressure shift coefiicient determined from the data. Line 
centers for 1^12 and U34 obtained in this way were then 
grouped as a function of sidereal time, where time zero 
has been set as the time in 1995 when we obtained our 
first data. 

Non-zero values for 5i'i2 and 5i'34 could arise from sys- 
tematic effects which lead to variations in the parameters 
affecting the line centers - particularly variations having 
a period of w 24 hr. Principal concerns are possible day- 
night variations of the magnetic field strength, and of the 
density and temperature of the Kr gas. 

Day-night variations of several °C in the temperature 
of the experimental hall lead to oscillations in the mag- 
netic field strength of the persistent-mode superconduct- 
ing solenoid of < 1 ppm. Changes of 0.05 ppm in the 
field strength were easily resolved, and the oscillation's 
effects on the line centers were accounted for in extracting 
the line centers. Temperature changes also affected the 
diamagnetic shielding constant of the water in the NMR 
probes used to monitor the field (the probes were not 
temperature-stabilized, but were in good thermal con- 
tact with the microwave cavity which was temperature- 
stabilized to 0.1 °C). The maximum conceivable 2 °C 
day-night changes in water temperature would change 
the NMR frequencies by 0.02 ppm. leading to errors in 
the line centers of about 2.5 Hz; of opposite sign for 1^12 
and 1^34. This potential effcK't is well below the statistical 
sensitivity for sidereal variations of 12 to 15 Hz. 

The effect on these data of the variation of Kr pressure 
with time has been evaluated. The front end window to 
the Kr stopping target was of 3 mil mylar, and flexed with 
day-night variations of the external atmospheric pressure. 
This induced fractional day-night oscillations in the Kr 
gas target pressure which were measured to be about 
2.5 X 10^''. Through pressure shift coefficients of about 
-16.5 kHz/atm for 1/12 and -19.5 kHz/atm for 1/34, the 
resulting shifts in the line centers (typically 7.5 Hz in 
1/34 and 6 Hz in P12) were automatically accounted for in 
performing the linear extrapolation to zero pressure, and 
should not contribute any significant time variation to 
Vij . The pressure shift coefficients depend on the average 
velocities of the atoms, and so are functions of tempera- 
ture. The fractional changes in the transition frequencies 
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with temperature (measured in hydrogen and its isotopes 
[6]) are roughly 1 x 10~^^ °C~^Torr~^. Given the temper- 
ature stability of the Kr gas of about 0.1 °C, temperature 
dependent errors in the extrapolation of the line centers 
to their vacuum values would be limited to a few Hz, well 
below the statistical sensitivity of our test. 

Other potential concerns involve the two frequency ref- 
erences used in the experiment - the proton precession 
frequency forming the basis of the magnetic field deter- 
mination, and the Loran-C 10 MHz frequency reference 
used for the NMR and microwave frequency synthesiz- 
ers. The Loran-C standard is based on hyperfine tran- 
sitions in Cs with mi?=0, and so is insensitive to any 
preferred spatial orientation, and would not introduce 
a signature for Lorentz violation into the spectroscopic 
measurements. Time dependent gravitational red or blue 
shifts of the frequency standard, caused by tidal distor- 
tions of the earth, are below the 10~^^ level and too 
small to be seen. Finally, bounds on clock comparisons of 
^^^Hg and ^'^^Cs [7,8] place limits on the Lorentz violat- 
ing energy shifts in the precession frequency of a proton 
of 10~^^ GeV, which imply the NMR measurements are 
free of shifts well below the Hz level. 

All the data obtained in 1995 and 1996 are plotted as 
a function of time measured as a function of a sidereal 
day in Fig. 3, where twelve points at « 2 hr. intervals 
are plotted, and the vertical scale is in Hz. The data for 
1/12 and 1/34 were fit separately by the functions 

i^ij (i) = {I'ij) + 5vij sin {2-Kt + (9) 

where t is the time in sidereal days, and the fit parameters 
arc {vij) ; the amplitude of the possible time variation 
5vij^ and the phase <pij (where no phase relation was 
assumed between y\2 and ^34). The amplitudes for 5ui2 
and (5i^34 are consistent with zero, — 4±13 Hz and — 13±15 
Hz respectively. 

As stated above, the theory being tested requires 
(5i^i2 ~ — <5i^34. A plot of z/12 — ^34 versus sidereal 
time is shown in Figure 4, and fit for a sinusoidal vari- 
ation (as in Eqn.9) where a common phase is assumed 
between P12 and z/34. The data exhibit no variation with 
time within ± 20 Hz, which corresponds to a 68% confi- 
dence level (one sigma) limit on the non-rotating frame 
components (see Eqn. 4) 

\J {b'^xY + (b'irf < 2 X 10-23 GeV. (10) 

The figure of merit of these results as a test of CPT 
violation is taken as 

^^m"+(Kfi-^, - « 5 X 10--. (11) 

The limits on Sv^^ and 15(1^34 — V12) yield similar values. 
Other choices for the denominator in the figure of merit 
such as Vij or Av are inappropriate since the former are B 



field dependent, and the latter is not a fundamental prop- 
erty of the muon. It is natural to take the muon mass 
as the basic energy scale since it is an input parameter 
of the standard model from which the atomic properties 
are derivable. The ratio of Lorentz and CPT violating 
deviations of muonium energy levels to the muon mass 
might be expected to be suppressed by the ratio of the 
low-energy (standard model) scale to the Planck scale 
of the underlying theory [9]. In this context the result 
achieved here, 5 x 10-^^, is comparable to the dimension- 
less scaling factor m,^jMp ~ . 

Bounds on a different linear combination of the muon 
parameters b^j,d^j^, and H'^j^ than from muonium can 
also be obtained from precision measurements of the 
anomalous magnetic moment of the muon = {g^—2) /2 
[5], where one can search for sidereal variations of the 
anomalous precession frequency (the difference be- 
tween spin and cyclotron frequencies). Also, the param- 
eter which is not tested in the latest muonium exper- 
iment, is bound at the level of 10^^^ GeV by the mea- 
surements of the anomalous precession frequencies of 
^+ and ^~ measured at CERN and BNL [5,10]. Further 
improvements in the Lorentz and CPT violating param- 
eters for the muon, coming from muonium, will require 
higher intensity muon sources, as the uncertainties are 
predominantly statistical. 

In conclusion, no unambiguous violation of CPT or 
Lorentz invariance has been observed although tests have 
been done on many systems since 1960 [11]. The results 
presented here represent Planck scale sensitivity, and are 
an order of magnitude improvement in sensitivity over 
previous limits on the muon parameters of the theory. 
They also represent the first test for sidereal variations 
in the 2nd generation of leptons, where Lorentz and CPT 
violating effects may be enhanced [5]. Higher sensitivity 
tests as well as tests on different physical systems will 
surely be made [3,12]. 

We thank V.A. Kostelecky for helpful discussions, and 
the U.S. DOE and BMBF (Germany) for supporting this 
research. 
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FIC. 1. Breit-Rabi energy level diagram of ground state 
muonium. At high fields, the indicated transitions, 1^12 and 
1/34 are essentially muon spin fiip transitions. 
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FIC. 2. Muonium resonance lines (data and fit) for 1^12 
taken using magnetic field sweep on the left, and microwave 
frequency sweep on the right, for muonium atoms which have 
decayed in selected time intervals after formation. 
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= <i^i2> + 51^,2 sin(27Tt+^i'2) 
= -4 +/- 13 Hz 
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2565762965 Hzlvs. Sidereal Time 




i^34(t) = <i'34> + (5l/34 Sin(27Tt + ^34) 

6i/,4 = -13 +/- 15 Hz 
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FIC. 3. Two years of data on 1^12 and F34 are shown binned 
versus sidereal time and fit for a possible sinusoidal variation. 
The amplitudes are consistent with zero. 
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- (1/34— i/ij) — 6682231 65 Hz vs. Sidereal Time 
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(f'34-l'l2)(t) = <l'34-J^12>+l5(l'34-i^12)sin(27Tt + ^) 

-6(1/34-1/,,) = -16 +/- 20 Hz 

-I iVr, , , ,i< < 
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A;/-4463302765 Hz vs. Sidereal Time 




Ai/(t) = <Ai/> + (5Ai/sin(27Tt + ^) 
SAu = -8 +/- 20 Hz 
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FIC. 4. Two years of data on 1^12-1^34 and i/i2 + 1^34 = Ai/ 
are shown binned versus sidereal time and fit for a possible 
sinusoidal variation. The amplitudes are consistent with zero. 
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